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ABSTRACT 

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned 
following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or 
more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities 
for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral 
energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) 
onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR 
wavelengths. We compare FIR-based SFRs to the more commonly used 24 yum and UV SFRs. We find that SFRs from 24 yum alone are higher by 
a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24 yum)> 12.2 Lq. The SFGs 
and AGNs tend to exhibit the same 24yum excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction 
correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and 
the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions. 

Key words. Galaxies:evolution - Galaxies :starburst - Galaxies:fundamental parameters - Cosmology:observations - Infrared: galaxies 
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1. Introduction 

The evolution of massive galaxies around z~2, when the cosmic 
star-formation rate density was at its peak, is currently the sub- 
ject of intense study. A key parameter in constraining their na- 
ture and stellar-mass buildup is the instantaneous star-formation 
rate (SFR). While the re are various methods for estimating SFRs 
(see iKennicutnl 19981 for a review), it is notoriously difficult to 
constrain SFRs of z~2 galaxies with an accuracy better than a 
factor of a few. One of the reasons is that calibrations between 
observed luminosities (rest-frame UV, mid-IR, and submm be- 
ing the most accessible and widely used) and bolometric SFRs 
are based on local galaxies and their validity at high z has not 
yet been firmly established. 

Both UV (BX/BM) and optical/near-IR (BzK) selected 
samples of star-forming galaxies at z~2 have been success- 
ful in producing large samples for galaxy evolution stud- 
ies ('Reddv et ak'^OOS*; 'Pere z-Gonzalez et al.l l2005l: iDaddi et al.1 
[2007; Santini et al. 2009). However, UV-based SFRs are very 
sensitive to dust obscuration and there are uncertainties about the 
applic able reddening law s and dust/source geometry at high z 
(e.g., iReddv et al.ll2010|) . Estimates from the 24 yum luminosi- 
ties rely on template-based extrapolations to the total IR lumi- 
nosity that are poorly constrained at high z and may be affected 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



by an AGN. Likewise, submm luminosities rely on assumptions 
about IR SEDs, which depend on poorly known dust tempera- 
tures, composition, and spatial distribution. 

Cosmological simulations predict a tight correlation between 
stellar mass and SF R with a slope similar to that observed (e.g., 
lFinlatoretani2006 '). but with SFRs that are a factor ~4 lower 
than those inferred for z~2 massive galaxies from a combination 
of Uy and 24 nm l uminosities dPaddi et al.ll2007HWilkins et aTl 
120081 : iDamen et al.1.200 9). Comparisons between the derivative 
of the mass density distribution and direct measurements of 
SFR density also lead to a discrepancy, in particular at z>2 
(e.g. Perez-Gonzalez et al. 2008). Suggestions on how to re- 
solve this discrepancy include changes to the adopted initial 
mass functions (IMF | Dave,i2008.) and star formation efficiency 
dKhochfar & Silkl l2009h . but these cannot be properly tested 
given the uncertainties associated with the common SFR indi- 
cators themselves. 

The photo-detect or array camera and spectrometer (PACS, 
Poglitsch et al. '2010^ onboard the Herschel space observatory 
(Pilbratt et al. 2010) offers a new opportunity to infer the SFR 
of z~2 galaxies from their far infrared (FIR) emission. PACS 
observes at 160 jim. (45-65 ;um rest frame for the sample used 
here), close to the peak of the emission from dust heated by 
young stars, away from AGN-heated dust emission and with un- 
precedented spatial resolution that reduces confusion noise. It 
allows the most reliable and least biased SFR measurements to 
date, which avoid large SED extrapolation errors and avoid the 
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uncertainties of attenuation corrections. In this letter, we report 
on the first measurements of FIR-based calorimetric SFRs of op- 
tically selected galaxies at 1.5<z<2.5 and compare them to re- 
sults from previous UV and 24 //m measurements. We adopt a 



Q„, = 0.3, Qa = 0.7, 
throughout this letter 



and Ho - 70 km s Mpc cosmology 



2. Data, sample and SFR(160jum) 

The PACS data consists entirely of PACS evolutionary probe 
(PEP) guaranteed-time observations in the GOODS-N field, 
taken during the science demonstration phase of the Herschel 
mission. For deta ils of the observat ions and data reduction pro- 
cess, we refer to iBerta et al.l (l2010h . Fluxes were extracted us- 
ing Spitzer-MWS 24 ^va sources with 3 cr detections or bet- 
ter as priors (Magnelli et al. 2009). At the current PACS depth 
(5.7 mJy, 3cr, at 160 /^m), testing showed that practically all 
PACS sources are detected at 24 //m. Hence, using 24 yum 
sources as priors allows a complete PACS source extraction. 
The priors ai-e matched to bviz (ACS), JHK (FLAMINGO), 
and Spitzer-IKAC p hotometry. Spectroscopic redshifts are from 
iBarger et al.l ( l2008l) a nd photometric redsh ifts were derived us- 
ing the code EAZY (JB rammer et alJl2008h . The Chandra 2Ms 
catalog ([Alexander et a l. 2003) is matched to the PACS sources 
and is used to flag X-ray sources. Using the iRanalli et alj 
(l2003h LJLfiR relation for SFG, the expected X-ray flux of a 
log(LIR/Lo)=12 galaxy at z=2 is almost an order of magnitude 
lower than the 1.4 x 10"'^ ergs cm"^ s"^ hard-band (observed 
2-8 keV) sensitivity limit of the 2 Ms catalog. We therefore 
treat every X-ray source with a hard-band detection as an AGN. 
One source without an X-ray detection but which exhibits a clear 
power-law SED in the IRAC bands, was also flagged as an AGN. 

In this study, we are interested in 1.5<z<2.5 range and re- 
quire redshift information. Hence, only sources for which we 
have spectroscopic redshifts, and/or were successfully fitted by 
EAZY for a photometric redshift are considered in the following 
selections. In addition, we require that in a 10" radius around a 
24 jum prior there is no more than one other 24 /vm source and 
it may have no more than 50% of the prior flux. This condition 
ensures that the sample does not experience significant flux con- 
fusion between MIR priors and/or FIR fluxes due to neighbors. 
These criteria result in a sample of 23 detected sources in total, 
17 classified as SFGs, and 4 of which have a spectroscopic red- 
shift. 306 24 /urn priors without a PACS detection are available 
for stacking analysis (see below). 

We calculate SFRs by assuming that they are proportional 
to the inte grated infrared luminosity (LIR, typically 8-1000 yum 
lKennicut l''l998'). We c onvert PACS 160 /um fl uxes into LIR by 
fitting die flux using the lCharv & Elbaj (1200 ll hereafter: CEOl) 
SED library, in which for a given redshift and band flux, a unique 
solution to LIR exists. The 100 //m filter is not used because 
at z~2, emission at rest-frame wavelengths may already be af- 
fected by hot dust and suffer the same systematic effects as the 
24 fim based SFR. In addition, at these redshifts 1 00 jum is on the 
same side of the SED peak and too close to 160 fiva to provide 
meaningful temperature constraints, without longer wavelength 
data. Combined PACS and SPIRE sub-mm data show that fitting 
CEOl SEDs to the 160 //m flux provides a reliable estimate of 
the total 8-1000 /v m luminosity (LIR) of the large-grain dust for 
z~2 galaxies (Elb az etal.ll2010l) . 

To explore the relationship between direct FIR luminosities, 
24 fim extrapolated luminosities, and UV-SFRs to lower lumi- 
nosities and SFRs, we also perform a stacking analysis. We stack 
the 160 yum image centered on the prior positions and measure 




Fig. 1. Total (8-lOOOA) IR luminosity from 160 //m versus that 
derived from 24 fim. Detections with spec-z are in blue, detec- 
tions with phot-z only are in red. SFGs are plotted as circles and 
AGNs as X marks. The error bars include only photometric er- 
rors. Squares represent mean luminosity of stacked SFG (green) 
and X-ray (black) undetected sources, with error bars indicat- 
ing the error in the mean luminosity. Horizontal bars beneath the 
stacks indicate the min-max values in the stack with the number 
of stacked sources above them. 



a mean flux for the stack. The mean luminosity of the stack, 
LIR(160 jum)5,flcA, is defined as the luminosity that will result 
in the 160 jum stacked flux, for the redshift distribution of the 
sources in the stack and the adopted CEOl SEDs. In this way, 
both k-corrections and luminosity distances are accounted for 
without bias. The error in the mean luminosity is calculated us- 
ing a bootstrap method. The same number of sources are resam- 
pled with replacements from the original stack to produce a new 
stacked image. The new mean luminosity is then calculated. The 
error is finally calculated from the distribution of luminosities in 
the repeated resampling. 



3. 24jum based SFR 

Spitzer-MIPS 24 yum fluxes are converted to LIR by fitting the 
24 //m flux to CEOl SEDs using the same method we use for the 
PACS 160 yum fluxes (Sect. E]). This SED library is very com- 
monly used in 24 fim (and other) studies where LIR is often 
derived from a single photometric point. Many other SED li- 
braries exist, but the vLv(8 /im)/LIR ratio that is relevant to this 
study is typically calibrated to match local galaxy observations 
and vary by factors smaller than those discussed below. At z~2, 
the rest-frame wavelengths are shorter than 10 //m and 24 //m 
fluxes probe the IR emission at the edge of the relevant range, 
far from the SED peak, making them sensitive to extrapolation 
errors. In addition, at these wavelengths polycyclic aromatic hy- 
drocarbon (PAH) emissions contribute significantly and the ra- 
tio of their fluxes to LIR may create significant scatter. PACS 
160 yum fluxes probe the emission much closer to the SED peak 
and are unaffected by the above issues. 

In Fig. [T] we plot LIR as derived from 160 jum versus that 
derived from 24 fj.m for 1.5<z<2.5. Individual sources are de- 
tected starting at LIR(160yL/m)> 6x10" Lq at redshifts close to 
z~1.5. The 24 yum sources without a PACS detection (and no X- 
ray counterpart) are stacked in luminosity bins as indicated by 
the bars in the figure. Stacked luminosities represent the mean 
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Fig. 2. The log ratio of 24 fim to 160 yum total IR luminosities 
(proportional to SFR) as a function of LIR from 24 yum. Colors 
and symbols are similar to Fig. [T] with black triangles represent- 
ing the mean of all (detections and non-detections) SFGs in the 
stacked range. The red arrow at top left indicates the typical un- 
certainties caused by phot-z errors. 



luminosity as defined in Sect.|2]and the stack error bars indicate 
the error in the mean LIR(160 yum) of the stack. Since we stack 
in luminosity and not flux bins, the typical redshifts of sources 
in a given stack are higher than those of detections with similar 
LIR(160 yum), which is one reason why not all sources in the lu- 
minosity bin are detected. It is clear that while there is some scat- 
ter, for most sources LIR(24 yum) is higher than LIR(160 yum). 

Figure 12] plots the LIR(24 yum) to LIR(160 /im) ratio as a 
function of LIR(24 yum). The detections appear to show a linear 
trend, but this is mostly due to a selection eff'ect. The combined 
mean for the detection and non-detections (detected in stacks) is 
much higher and the individual detections represent a fraction of 
the distribution for which at a given LIR(24 yum), LIR(160 fim) 
is high enough to be detected. There is a weak 24 yum lumi- 
nosity dependence and the mean of the excess increases with 
LIR(24 yum), at least for log LIR(24 yum)> 12.2. Below this lu- 
minosity, it is harder to determine the slope, but the stack indi- 
cates that it may flatten out. The dashed and dotted lines in Fig.|2] 
are ad hoc fits to the binned combined-means slope above and a 
constant off'set below LIR(24)=12.2 



0.38 logLIR(24)- 4.0 
0.6 



logLIR(24)> 12.2 
logLIR(24)< 12.2 



(1) 



Though a large fraction of the sample relies on photometric 
redshifts, errors in redshifts will move the points in LIR(24yum), 
but will have a much smaller effect on the LIR ratio. The red ar- 
rows at the top-left of Fig. |2] demonstrate how the typical phot-z 
errors, propagated into the quantities plotted in the figure, will 
affect the values. The photo-z induced dispersion is directed 
along the diagonal trend of detected sources and probably con- 
tributes to this pattern with the selection effects. 

X-ray AGNs do not appear to differ from the rest of the sam- 
ple in Figs. [Hand |2] except perhaps for the one source with the 
highest luminosity. This source also exhibits a very clear power- 
law SED in the mid-IR. For sources completely dominated by 
the AGN, fitting a CEOl SED to the 24 /im emission is clearly 
inappropriate. Even though we can only produce one AGN stack 
that results in a good detection (Fig. [TJ, the mean luminosity 
is similar to that of the SFG stacks. This can be interpreted in 
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Fig. 3. The sum of SFR from 160 jum and Wobs (no attenuation 
correction) versus SFR derived from attenuation-corrected UV 
flux. Blue symbols are for spec-z, red are for phot-z. Circles are 
SFG detections, x-marks are AGNs. Errors in detected sources 
are photometric only. Squares are mean SFR for stacks with the 
error in the mean. Horizontal bars below indicate the min-max 
range of values in the stack with the number of stacked sources 
noted above. 



two ways: either the flux at 6.8 < A < 9.6 yum rest frame is 
dominated by the starburst emission at all the luminosities con- 
sidered here, or the flux excess in all sources is due to an AGN, 
which is obscured in most. The latter claim and the relatively 
tight LIR(24yum)/LIR(160yum) ratio would imply a tight relation 
between the AGN luminosity and the galaxy 's starburst com- 
ponent FIR luminosity, which is not observed dLutz et al.ll2010l : 
Shao et al. 2010). We conclude that while a hidden AGN may 
contribute to the scatter, it is not likely to be the main cause of 
the general 24 fj.m excess. Enhanced emission (relative to local 
galaxy SEDs) in PAH features between 6.2-8.6 //m, that enter 
the 24 yum filter at z= 1 .5 is a more plausible explanation. 

These new resu lts corroborate the main findings of 
iMurphv et all (l2009h . based on Spitzer 16, 24, 70 yum, and 
SCUBA 850 fim observations. For z~2 galaxies, they find that 
SFRs using only 24 yum are overestimated by a factor of 5 and 
conclude t hat a hidden AG N can only account for a fraction of 
this excess. IPapovich et al.l 12007), using Spitzer 24, 70 and 160 
yum stacked fluxes also find that LIR(24 yum) is in excess. At the 
high end, they find 24 yum overestimates LIR by a factor of 2- 
10. They too conclude that X-ray sources have a 24 yum excess 
similar to SFGs without X-ray counterparts. 



4. UV-based SFR 

The UV continuum at wavelengths longer than the Lyman edge 
can be used as a SFR indicator. The UV emission is dom- 
inated by young, massive, and short-lived stars and a SFR 
that is required to sustain their population can be calculated. 
The SFR(UV) is sensitive to the attenuation correction used, 
which varies significantly in the UV between different extinction 
curves, all derived from local galaxies. Testing the accuracy of 
SFR(UV) is essentially a test of the attenuation correction. Here 
we use B and z band photometry (probing rest frame ~1500 A 
and ~3500 A at z~2) to estimate their SFR and dust attenua- 
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Fig. 4. The log ratio of UV SFR to the combined LIR and \5Yohs 
SFR versus Ks magnitude. Colors and symbols are the same as 
in Fig. [3] 



tio n. For the optical e xtinction, we apply the calibration given 
bv lDaddiet all (12004 



E{B - y) = 0.25(B -z + 0.1)ab. 



(2) 



The effective attenuation at 1500 A is obtained using the 
ICalzetti et all (12000.) extinction law A1500 = \QE{B - V), which 
was derived for local starburst galaxies. 

Th e BzK sample is selected according to the iDaddi et al] 
( |2004|) color criterion. We require at least 3 cr detection in all 3 
bands, a cut at Ks<22 AB, and a redshift (spectroscopic or pho- 
tometric) 1.5<z<2.5. PACS 160 yum fluxes were extracted using 
24 /zm priors, which implies a 24 jjm. detection for this selection. 
It also means that the BzK sample is a subsample of that used 
in Sect. |3] The sample includes 15 sources detected at 160 /^m, 
11 classified as SFG. Seventy-five sources without 160 /^m de- 
tections are used for stacking. Strictly speaking, the requirement 
for a 24 ^im prior is irrelevant to both SFR(UV) and SFR(160) 
and may even bias the selection against sources with lower ob- 
scuration. To improve the stacked signal, we also include BzK 
sources without a 24 /^m detection, but verify that introducing 
this requirement gives consistent results. In converting from lu- 
minosities to SFRs, we follow the Kennicutt ( 1998) relations, but 
lower the SFR by a factor 1.7 to convert from a Salpeter initial 
mass function (IMF) to that of .Chabrier (.2003.) . though for the 
results discussed here the choice of IMF makes little diff'erence: 



SFRuAMJyr-'] = 
SFRuvW^lyr-'] = 



LIR 



9.85 X 109 Lo' 
Lv(1500A) 



1.5 X 1028 ergs- iHz- 



(3) 
(4) 



In Fig. |3] we plot the combined SFR from LIR(160 jim) and the 
observed (uncorrected) UV flux versus the attenuation-corrected 
UV SFR, SFR(UVcorr)- The LIR is assumed to be caused by the 
reprocessing of radiation from young stars by the dust and the 
addition of the SFR from the observed UV luminosity accounts 
for the escaped part of the UV radiation. In our sample, the lat- 
ter makes only a small contribution with a mean of 2.4% and 
maximum of 8% of the total SFR for PACS sources. Our de- 
tection limit with PACS is at about 70 Mq yr', below which 
we can only study sources by means of stacking. However, for 
a cut level Ks < 22 the small number of BzK galaxies below 
the PACS threshold limits the stacking possibilities. The de- 
tections exhibit good agreement between SFR(UVco/-/-) and the 



combined SFR(LIR-i-UVofc5). The mean log-ratio of the detec- 
tions is -0.02 dex and standard deviation of 0.25 dex. We stack 
PACS non-detections by SFR(UV„,rr) bins. For the low SFR 
stack, the two SFR estimators are in very good agreement, while 
for the other, the value of SFR(UV corr) is slightly too high. The 
overall mean log-ratio of stacks and detections is 0.3 dex and the 
standard deviation is 0.35 dex. We therefore conclude that when 
using the Calzetti UV extinction law for SFR(UV)>40 Mq yr"' 
the SFR is overestimated by a mean factor ~2 with a scatter of a 
similar magnitude. 

The scatter in Fig. [3] exhibits a systematic pattern. For 
SFR(UVc„rr)>100 Mq yr-1, this estimator tends to overpredict 
the SFR of the PACS detections by a factor~1.5 and below this 
rate to underpredictby roughly the same factor. This trend is also 
seen in Fig|4] where the SFR indicator ratio is plotted vs. the Ks 
magnitude. Since the pattern appears to be also a function of Ks 
magnitude (and of both B and z mag to a lesser degree), it is 
unlikely to be caused by phot-z errors. It is not preserved when 
plotting against either redshift or extinction. Both stacks for the 
Ks magnitude bins lie above the zero-line, indicating a spread 
in the full population and shift the overall mean upwards. The 
difference between stacks and detections is likely a selection ef- 
fect: for a given SFR(UV) (or Ks), one will generally tend to 
detect in 160 yum galaxies for which the attenuation is under- 
predicted. However, the clear dependency on both SFR and Ks 
magnitude and the observed trend for the stacks might suggest 
there is more to it. The geometry an d metallicity of the galaxies 
may play a role in this. IReddv et al.l (12010.) suggested that some 
young galaxies may have a different UV reddening law, although 
these galaxies represent only 13% of their sample. The current 
sample is too small to draw conclusions about these effects. 

Four X-ray AGNs that are selected by BzK were detected 
with PACS at 160 /im. Three are indistinguishable from the SFG 
in Figs. |3] & |4] and both their UV and FIR emission must be 
dominated by the star formation. The fourth (log-ratio=-0.75) is 
an extremely bright AGN and both SFR indicators are probably 
not applicable to this case. This is the only 160 /im source that 
does not appear in Figs. [1] and |2] as it is too bright in 24 yum 
to be fit by CEOl SEDs. There were not enough non-detected 
AGN-BzKs to produce a meaningful stack. The presence of an 
AGN, except in the extreme cases, does not significantly affect 
the SFR(UV). 
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